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Abstract 

Background: Resequencing can be used to identify genome variations underpinning many morphological and 
physiological phenotypes. Legume model plant Medicago truncatula ecotypes Jemalong A17 (J. A17) and R108 
differ in their responses to mineral toxicity of aluminum and sodium, and mineral deficiency of iron in growth 
medium. The difference may result from their genome variations, but no experimental evidence supports this 
hypothesis. 

Results: A total of 12,750 structure variations, 135,045 short insertions/deletions and 764,154 single nucleotide 
polymorphisms were identified by resequencing the genome of R108. The suppressed expression of MtAACT that 
encodes a putative aluminum-induced citrate efflux transporter by deletion of partial sequence of the second intron 
may account for the less aluminum-induced citrate exudation and greater accumulation of aluminum in roots of 
R108 than in roots of J. A17, thus rendering R108 more sensitive to aluminum toxicity. The higher expression-level 
of MtZpt2-l encoding a TFIIIA-related transcription factor in J. A17 than R108 under conditions of salt stress can be 
explained by the greater number of stress-responsive elements in its promoter sequence, thus conferring J. A17 
more tolerant to salt stress than R108 plants by activating the expression of downstream stress-responsive genes. 
YSLs (Yellow Stripe-Likes) are involved in long-distance transport of iron in plants. We found that an YSL gene was 
deleted in the genome of R108 plants, thus rendering R108 less tolerance to iron deficiency than J. A17 plants. 

Conclusions: The deletion or change in several genes may account for the different responses of M. truncatula 
ecotypes J. A17 and R108 to mineral toxicity of aluminum and sodium as well as iron deficiency. Uncovering 
genome variations by resequencing is an effective method to identify different traits between species/ecotypes that 
are genetically related. These findings demonstrate that analyses of genome variations by resequencing can shed 
important light on differences in responses of M. truncatula ecotypes to abiotic stress in general and mineral stress 
in particular. 

Keywords: Resequencing, Medicago truncatula, Aluminum toxicity, Aluminum- activated citrate transporter, Salt 
stress, MtZpt2-l, Iron deficiency, Yellow Stripe-Likes 



Background 

Legume is the second most important crop family in the 
world, and is one of primary sources for the consump- 
tion of human and animals [1,2]. Acquisition of nut- 
rients from soil is a prerequisite for plant growth and 
development. Plants are frequently exposed to adverse 
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mineral stress in soils, including aluminum toxicity in 
acid soil, salt stress in saline soil and iron deficiency in 
alkaline soil. Plants have evolved numerous mechanims 
to adapt to these stressed environments [3-5]. Under- 
standing of the molecular mechanims by which plants 
respond and adapt to the mineral toxicity and deficiency 
is a major challenge in modern plant biology. 

As a model legume species, Medicago truncatula 
Gaertn has been widely used to study functional genomics 
because of its small diploid genome, self-fertility, short 
generation cycle and easy transformation [6]. There are a 
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number of ecotypes of M. truncatula with large genetic 
variations [7]. Of the ecotypes, M. truncatula ecotype 
Jemalong A17 (J. A17) has been used for the whole- 
genome sequencing and physiological studies [8-10], 
while ecotype R108 is often used for gene transform- 
ation because of its superior in vitro regeneration [11]. 
M. truncatula ecotype R108 differs from its counterpart 
J. A17 in traits associated with development, and biotic/ 
abiotic responses. For example, treatments of J. A17 with 
methyl jasmonate and ethylene induce resistance to fungal 
pathogen Macrophomina phaseolina, while these treat- 
ments fail to induce resistance in R108 to the fungal 
pathogen [12]. In addition, rhizobial-induced expression 
of chitinase gene between the two ecotypes is also dif- 
ferent [13]. The two ecotypes exhibit different tolerance 
to salt stress, such that ecotype J. A17 is more tolerant 
to salt stress than R108. Further studies reveal that a 
TFIIIA-related transcription factor gene, MtZpt2-l 
shows different expression in the two ecotypes, and that 
overexpression of MtZpt2-l in roots confers enhanced 
tolerance to salt stress [14,15]. Our previous work re- 
vealed that the two ecotypes also differed in their toler- 
ance to deficiency in mineral nutrients. For example, 
ecotype J. R108 was more sensitive to iron deficiency 
than ecotype J. A17 [16]. Despite the morphological and 
physiological differences between the two ecotypes, few 
studies have investigated the molecular mechanisms 
underlying the differences due to lack of information on 
the genome of R108. 

DNA sequences contain all the genetic information, 
and genome variations such as structure variations (SVs), 
short insertions/deletions (indels) and single nucleotide 
polymorphisms (SNP) can explain many variations in 
morphological, physiological, and ecological traits [17-21]. 
Resequencing technology provides a powerful tool to 
study these variations among species/ecotypes that are 
closely related genetically. For instance, Thellungiella sal- 
suginea exhibits exceptionally high resistance to cold, 
drought, and oxidative stresses as well as salinity [22-24]. 
The number of members in gene families with known 
functions associated with responses to abiotic stresses in 
T. salsuginea is greater than in Arabidopsis thaliana, in- 
cluding those gene families of RAV, NF-X1, GRAS, HSF, 
HKT, CIPK and CDPK [25]. Furthermore, it has been re- 
ported that maize inbred-line Mol7 exhibits eminent 
heterosis due to its deletion of eighteen genes [19]. The 
two widely used M. truncatula ecotypes Jemalong A17 
(J. A17) and R108 have been reported to differ in their 
tolerance to salt stress [14,15] and iron deficiency [16]. 
To test whether the genome variations between J. A17 
and R108 may account for the differences in their re- 
sponses to mineral toxicity of aluminum and sodium 
and mineral deficiency of iron in growth medium, genome 
variations of M. truncatula ecotype R108 were analyzed 



by mapping the reads obtained from resequencing of 
R108 to the reference genome of ecotype J. A17. 

Results 

Response of J. A17 and R108 to Al 3+ and Na + toxicity, and 
Fe deficiency 

To examine the effect of Al 3+ on root elongation of the 
two ecotypes, the relative root elongation was determined. 
As shown in Figure la, root elongation was inhibited upon 
exposure of the two ecotypes to solution containing Al 3+ , 
and the Al 3+ -induced inhibition of root elongation in 
R108 was greater than in J. A17. Moreover, Al contents in 
R108 roots were higher than in J. A17 roots (Figure lb), 
implying that an exclusion mechanism may operate in 
ecotype J. A17 plants. Exudation of organic anions includ- 
ing malate and citrate to complex toxic Al 3+ in the rhizo- 
sphere is an important mechanism to tolerate Al [3,26,27] . 
Therefore, we monitored exudate of malate and citrate 
from roots of the two ecotypes in response to Al 3+ treat- 
ment. There was an increase in citrate exudation from 
roots of J. A17 and R108 plants by exposure to Al 3+ , and 
the Al 3+ -induced citrate exudation from roots of J. A17 
was greater than that of R108 plants (Figure lc). In con- 
trast to citrate, no significant increases in malate exud- 
ation from roots of the two ecotypes by exposure to Al 3+ 
were detected (data not shown). These results suggest that 
higher exudation of citrate may underpin the greater toler- 
ance of J. A17 to Al than R108 plants, de Lorenzo et al 
found that R108 is more sensitive to salt stress than J. A 17 
plants as evidenced by less suppression of root growth in 
J. A17 plants than in R108 plants [14]. In addition to root 
growth, Na + /K + ratio is an important indicator for toler- 
ance of plants to salt stress. Excessive accumulation of 
toxic Na + in plant cells, particularly in the cytosol, disrupts 
K + homeostasis, leading to dysfunction of plant cells, thus 
plants displaying high tolerance to salt stress often mi- 
nimize Na + uptake and/or maximize K + acquisition to 
maintain a low Na + /K + ratio [28]. Therefore, we compared 
the effect of salt stress on Na + and K + concentrations in 
the two ecotypes. No differences in both Na + and K + con- 
centrations in shoots of the two ecotypes were found 
when they were grown in the control medium (Figure 2a 
and b). When they were exposed to solution containing 
NaCl, an enhanced accumulation of Na + in both ecotypes 
was observed (Figure 2a). However, exposure to salt stress 
led to reductions in K + concentrations in shoots of both 
ecotypes, and the salt stress-induced reduction in K + 
concentration was greater in R108 than in J. A17 plants 
(Figure 2b). This led to an increase in Na + /K + ratio in both 
ecotypes, and the increase was significantly less in J. A 17 
plants than in R108 plants (Figure 2c). 

Our previous work showed that the ecotype J. A17 
was more tolerant to Fe deficiency than R108 by effi- 
ciently mobilizing Fe in the rhizosphere and transporting 
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Figure 1 (See legend on next page.) 





-Al +AI 

Treatments 



Wang et al. BMC Plant Biology 2014, 14:122 
http://www.biomedcentral.com/1471-2229/14/122 



Page 4 of 1 1 



(See figure on previous page.) 

Figure 1 Effect of Al 3+ on root elongation, citrate exudation and Al content in roots of J. A17 and R108 plants. The relative root 
elongation was determined by exposing 3-d-old seedlings of J. 17 and R108 to 5 uM AICI 3 (pH 4.5) for 2 days (a). Data are mean±s.e. with 
n = 10. Al contents in roots of J. A17 and R108 plants before and after exposure to 5 uM AICI 3 (pH 4.5) for 2 days (b). Data are mean±s.e. with 
n = 4. Citrate exudation rate from roots of J. A1 7 and R1 08 plants treated with 5 uM AICI 3 (pH 4.5) for 24 h (c). Data are mean ± s.e. with n = 5. 
* and ** indicate significant difference between genotypes within a given growth condition (-AI or + Al) at P< 0.05 and P< 0.01, respectively. 



of Fe from roots to shoots in J. A17 plants [16]. A simi- 
lar result showing that ecotype J. A17 had higher foliar 
Fe contents than R108 when grown in Fe-deficient me- 
dium was observed in the present study (Figure 2d). 
These results show that the two ecotypes differ in their 
tolerance to toxicities of Al and Na as well as Fe defi- 
ciency by differently regulating citrate exudation, Na up- 
take and Fe transport, respectively. 

Resequencing of R108 

Paired-end sequencing method was employed to rese- 
quence the genome of M. truncatula ecotype R108, and 
about 4.64 Gb original sequencing data were generated. 



High-quality reads of 4.28 Gb were obtained after ini- 
tially processing. The genome of R108 is 17% smaller 
than that of J. A17 [29]. This led to a sequencing mean 
coverage and depth of approx. 72% and 11 -fold over the 
whole genome, respectively (Figure 3 and Additional file 1: 
Figure SI). The coverage of chromosome 5 was the 
greatest among the chromosomes. In addition, we found 
a similar coverage of chromosome 5 in a tetraploid Medi- 
cago falcata (unpublished results), suggesting that Chr 5 
may be the most conserved in the genus of Medicago. 

Structure variations (SVs), short insertions/deletions 
(indels) and single nucleotide polymorphisms (SNPs) were 
identified by aligning the high-quality sequences against 




Treatments Treatments 
Figure 2 Effects of salt stress and iron deficiency on Na + and K + concentrations, Na + /K + ratio, and Fe concentrations in shoots of 
J. A17 and R108 plants. Concentrations of Na and K and Na + /K + ratio in shoots treated with and without 100 mM NaCI for 5 days were shown 
in panel (a), (b) and (c), respectively. Data are meant s.e. with n = 4. Fe concentration in shoots of 5-d-old seedlings of J. A17 and R108 plants 
exposed to control, Fe-sufficient medium (100 uM Fe-EDTA, +Fe) and Fe-deficient medium (1 uM Fe-EDTA) for 5 days (d). * and ** indicate 
significant difference between genotypes within a given growth condition at P< 0.05 and P< 0.01, respectively. 
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Figure 3 The sequencing coverage of 8 chromosomes in the genome of R108 against to the reference of J. A17 genome. One hundred 
kb was defined as one window. 



the reference genome of J. A17. We obtained a total of 
12,750 SVs, 135,045 indels and 764,154 SNPs in the gen- 
ome ofR108 (Table 1). 

Structure variations are important types of differences 
among individuals of the same species, and can cause 
large alterations to the genome, resulting in the dif- 
ferences in phenotypes. We identified 10,964 deletions, 

Table 1 The number of SVs, indels and SNPs in the R108 
genome 





Genome variations 


Numbers 


sv 


Insertion 


1,239 




Deletion 


10,964 




Others 


547 


Indel 


Insertion 


67,087 




Deletion 


67,958 


SNP 


Homozygosity 


660,168 




Heterozygosity 


103,986 



The data were obtained by mapping the reads of R108 obtained from 
resequencing to the reference genome of ecotype J. A17. SVs were separated 
into insertions, deletions and others. Indels were separated into insertions and 
deletions. SNPs were separated into homozygosity and heterozygosity. 



1,239 insertions and 547 other SVs such as duplication, 
inversion and transposition by resequencing (Table 1). 
The quantity of deletions was more abundant than other 
SVs. This result is consistent with the forecast as the 
genome of R108 has been reported to be smaller than 
that of J. A17 [29]. We also identified 135,045 indels 
ranging from 1-5 bp in length. Among these short 
indels, the number of insertions and deletions was al- 
most equal (Table 1). Insertion of one bp and deletion of 
one bp were the mostly observed insertions and dele- 
tions, respectively, accounting for more than half of the 
total number of insertions and deletions (Additional 
file 1: Figure S2). Generally, genome variations were 
mainly accounted for by SNPs. Eighty-six percent of 
SNPs were homozygous over the whole genome (Table 1). 
For the SNPs within coding sequences, there were 70,695 
nonsynonymous and 57,124 synonymous SNPs, respect- 
ively. This led to a ratio of nonsynonymous to synony- 
mous nucleotide (Nonsyn/Syn) of 1.24. A similar ratio has 
been reported in soybean and rice, while the ratio in 
Arabidopsis (0.83) is smaller than our finding in the 
present study [17,18,30]. 
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Variations of mineral element-related genes 

The resequencing data obtained from M. truncatula 
ecotype R108 revealed that some genes involved in ac- 
quisition of mineral elements were deleted in R108 plants 
compared to ecotype J. A17 (Table 2). 

Aluminum-induced exudation of citrate from roots 
that is mediated by membrane transporters can detoxify 
toxic Al 3+ in the rhizosphere by forming non-toxic Al- 
citrate complex [3,26,27]. Several genes encoding the 
transporters of Al-induced citrate exudation belonging 
to MATE family have been identified [31-34]. Our rese- 
quencing data show that 771 bp in the second intron of 
a gene encoding a putative aluminum-activated citrate 
transporter (MtAACT) was deleted (Table 2, Additional 
file 1: Figure S3). The sequence of MtAACT was si- 
milar with the known Al-induced citrate transporters 
(Figure 4a). Expression-level of MtAACT in J. A17 was 
higher than in R108 in the absence of Al 3+ , and it was 
up-regulated in both ecotypes by exposure to A1C1 3 
with the Al-induced expression of MtAACT in R108 
being lower than in J. A17 plants (Figure 4b). 

Previous studies have shown that R108 is more sen- 
sitive to salt stress than J. A 17, and real-time qPCR 
showed that expression of MtZpt2-l is greater in J. A17 
than in R108 plants [14]. Overexpression of MtZpt2-l in 
roots of the salt-sensitive ecotype of M. truncatula con- 
fers enhanced tolerance to salt stress, suggesting that 
differential expression of MtZpt2-l is responsible for the 
difference in adaptation to salt stress. Resequencing al- 
lowed us to analyze the promoter sequence of MtZpt2-l 
in R108 plants. Stress-responsive-related ds-elements 
were identified by PLACE database between J. A17 and 
R108 (Figure 5). The number of MYC and W-box ele- 
ments was greater in J. A17 than in R108, which may 
underpin the higher expression levels of MtZpt2-l in J. 
A17 than in R108 plants under conditions of salt stress. 

YSLs (Yellow Stripe-Likes) are involved in long- 
distance transport of Fe in plants [35,36]. There are five 
YSLs in the genome of M. truncatula according to 
Mt3.5 assembly of the reference genome. Our resequen- 
cing results show that an YSL gene (Medtrlg007540) 
was deleted in the genome of R108 (Table 2). The pro- 
tein encoded by Medtrlg007540 is highly similar to Ara- 
bidopsis AtYSL3 (At5g53550) (Figure 6). The deletion of 
the YSL gene in the genome of R108 may account for 
the less accumulation of iron in the shoots of R108 
(Figure 2d). 



Discussion 

Identification of genome variations using resequencing 

Analyses of gene expression by methods such as transcrip- 
tome, microarray and DGE have been used to decipher 
the differential responses among species and cultivars/ 
ecotypes with close genetic background to abiotic stresses 
[14,37-40]. However, these methods are less effective when 
several mechanisms underlie the different responses to 
abiotic stresses. Moreover these methods cannot be used 
to analyze ds-acting regulatory elements. In contrast, 
resequencing technology can identify genome varia- 
tions which are responsible for morphological and physio- 
logical differences [41]. In addition, ds-acting regulatory 
sequences obtained from the resequencing can be used to 
pinpoint the differential expression in response to abiotic 
stresses. Estimation of phylogenetic relationships among 
Medicago species by genome resequencing has been re- 
ported [42] . However, genome resequencing has not been 
used to investigate responses of Medicago species to abi- 
otic stresses in general and mineral stresses in particular 
so far. In the present study, we utilized this technology to 
decipher the mechanisms underlying the different re- 
sponses of two M. truncatula ecotypes to aluminum tox- 
icity, salt stress and iron deficiency. 

Tolerance of M. truncatula to Al is achieved by citrate 
exudation 

Aluminum is the most abundant metal in the earths 
crust. Phytotoxic Al 3+ is solubilized when soil becomes 
acidified. Inhibition of root elongation is one of the ear- 
liest and most distinct symptoms exhibited by plants 
suffering from Al toxicity [43]. Plants have evolved nu- 
merous mechanisms to adapt to Al toxicity. Exudation 
of organic anions from root apices to chelate toxic Al 3+ 
in the rhizosphere is an effective way to detoxify Al tox- 
icity, thus conferring tolerance to Al toxicity [3,26,27]. 
Several Al-activated citrate transporters have been shown 
to be involved in regulation of Al tolerance. For instance, 
SbMATE in sorghum (Sorghum bicolor) and HvAACTl 
in barley (Hordeum vulgare) that belong to the multi- 
drug and toxic compound exudation (MATE) family 
have been identified to mediate Al-activated citrate exud- 
ation. Heterologous expression of SbMATE in Arabidopsis 
and HvAACTl in tobacco leads to enhanced citrate efflux, 
thus conferring tolerance to Al toxicity [31,32]. The ho- 
mologs in Arabidopsis and maize have subsequently been 
cloned [33,34]. 



Table 2 Deleted genes related to mineral stress in the genome of M. truncatula ecotype R108 plants 



Genes 


Positions in reference 


Annotations 


Positions of deletions 


Medtr8g036660 


MtChr8: 8363576-8369153 (-) 


Putative aluminum activated citrate transporter 


MtChr8: 8367361-8368131 


Medtr1g007540 


MtChrl: 996073-999271 (-) 


Metal-nicotianamine transporter YSL3 


MtChrl :980798-1 01 0225 



The data were obtained by mapping the reads of R108 obtained from resequencing to the reference genome of ecotype J. A17. 
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Figure 4 Similarity of MtAACT protein to other known AACT proteins and effect of Al 3+ on expression of MtAACT in J. A17 and 
R108 plants. Phylogenetic tree of known and putative Al-activated citrate transporters was constructed by MEGA 5 in panel (a). The accession 
numbers of SbMATE, HvAACTI, AtMATE, ZmMATEl, MtMATE, At4g38380 and At2g38330 in GenBank are ABS89149.1, BAF75822.1, NP_974000.1, 
ACM47311.1, XP_003627698.1, NPJ95551.5 and NP_181367.2, respectively. The expression of MtAACT in roots of J. A17 and R108 plants under 
the conditions of with or without 5 uM A1CI 3 (pH 4.5) in medium for 1 days (b). Data are mean±s.e. with three biological replicates. * and ** 
indicate significant difference between genotypes within a given growth condition (-AI or + Al) at P< 0.05 and P < 0.01, respectively. 



In the present study, we found that root elongation of 
R108 was more inhibited by Al than that of J. A17 plants 
(Figure la), suggesting that R108 is more sensitive to Al 
than J. A17. We uncovered deletion of partial sequence 
in the second intron of a gene encoding a putative Al- 
activated citrate transporter (MtAACT) in R108 plants 
by resequencing (Table 2). The amino acid sequence of 
this transporter is similar with the known Al-activated 
citrate transporters in other plant species (Figure 4a). In 
addition, expression of this gene was up-regulated by Al 
in both ecotypes with the magnitude of Al-induced ex- 
pression of MtAACT in R108 less than in J. A17 plants 
(Figure 4b), suggesting that expression of MtAACT is 
sensitive to Al 3+ . The suppressed expression of MtAACT 
in R108 relative to that in J. A17 plants is likely to be 
accounted for by the deletion of the partial sequence of 
the second intron. The intron deleted in R108 plants 
may activate gene expression by enhancers within in- 
trons and/or regulating chromatin remodeling. Several 
introns in plants are reported to increase the expression 
of genes. For example, the second intron from Arabi- 
dopsis agamous gene can function in both orientations 
to drive expression of reporter gene from a minimal pro- 
moter [44]. The first intron of Arabidopsis gene encod- 
ing elongation factor eEF-lp has the similar function to 



enhance gene expression [45]. The lower abundance of 
MtAACT transcripts in R108 than J. A17 when exposed 
to solution containing toxic Al 3+ may explain the less 
citrate released from roots of R108 than J. A17 plants in 
response to Al treatment (Figure lc). The reduced cit- 
rate exudation from roots of R 108 plants due to reduced 
expression of MtAACT would render R108 plants less 
effective to complex toxic Al 3+ in the rhizosphere, thus 
making it less tolerant to Al than J. A17 plants. The 
greater accumulation of Al in roots of R108 than in 
those of J. A17 is in line with this argument (Figure lb). 

Promoter analysis of MtZpt2-1 

The ecotype J. A17 plants have been shown to be more 
tolerant to salt stress than R108 plants [14,15]. A gene 
encoding a TFIIIA-related transcription factor, MtZpt2-l 
has been identified by its greater up-regulation in J. A17 
than R108 plants under salt stress [14]. MtZpt2-l can 
active the expression of many stress-responsive genes 
[46]. Several stress-related ds-elements were found by 
analyzing the promoter sequences of MtZpt2-l in both 
ecotypes (Figure 5). These stress -related ds-elements in 
MtZpt2-l can allow this gene to be up-regulated in re- 
sponse to abiotic stresses, thus participating in the regu- 
lation of tolerance to abiotic stresses. Two MYB-core 
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Figure 5 Analysis of MtZpt2-l promoter sequence of J. A17 and R108. The arrows above line represent c/'s-elements of J. Al 7, and that of 
below the line indicate elements of R108. 
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Figure 6 Sequence analysis of YSL protein family. Phylogenetic tree of these proteins was constructed by MEGA 5. The corresponding IDs 
were shown in the figure. 



elements and ABA-responsive elements (ABRE) have been 
shown to be involved in responses to osmotic stress and 
ABA, respectively [47,48]. However, the two ecotypes dif- 
fered in their promoter sequences of MtZpt2-l, such that 
MtZpt2-l of J. A17 plants had one more ACGT element, 
four more MYC elements and two more W-box elements 
than R108 plants. There are reports showing the involve- 
ments of these ds-elements in stress response [49-51]. 
The greater number of ds-elements of MtZpt2-l in J. A17 
plants may explain higher expression of MtZpt2-l in J. 
A17 plants than in R108 plants under conditions of salt 
stress, thus conferring their tolerance to salt stress. 

Function of YSL in iron transport 

YS1 (Yellow Stripe 1) has been identified to be involved 
in uptake of iron from soil by roots in maize [52,53]. 
Based on their sequence similarity to the maize YS1, 
eight YSLs (Yellow Stripe-Likes) were identified in Ara- 
bidopsis. AtYSLl, AtYSL2 and AtYSL3 are expressed 
most strongly in the vascular parenchyma cells [36,54]. 
The ysllysl3 double mutant displays strong interveinal 
chlorosis, and has reduced foliar iron content [36]. These 
findings suggest that YSLs act as key mediators in unload- 
ing iron to mesophyll cell after iron is transported from 
roots through xylem in plants [55]. 

Five YSLs were identified in the genome of Medicago 
according to Mt3.5. However, in the genome of R108, an 
YSL gene (Medtrlg007540) was deleted (Table 2). The 
protein encoded by the YSL gene had high similarity 
with AtYSL3 of Arabidopsis (Figure 6). We hypothesize 
that this protein may be involved in unloading of iron 
from the vascular tissues to mesophyll cells. The dele- 
tion of this gene in R108 plants would impair iron 
unloading to mesophyll cells, thus leading to the re- 
duced iron contents in shoots of R108 plants when 
grown in iron-deficient medium (Figure 2d). 

Conclusions 

The two M. truncatula ecotypes Jemalong A17 and 
R108 differed in their sensitivity to aluminum toxicity, 



salt stress and iron deficiency. Resequencing of M. trun- 
catula ecotype R108 uncovered a total of 12,750 SVs, 
135,045 indels and 764,154 SNPs by comparing with the 
reference genome of J. A17. We found that the partial 
sequence of the second intron of MtAACT that encodes 
a putative Al-activated citrate transporter was deleted. 
This partial deletion may lead to the lower expression 
level of MtACCT in R108 plants than that in J. A17 
plants in the absence and presence of toxic Al in the 
growth medium. The reduced expression of MtAACT in 
R108 plants in turn may render less exudation of citrate 
form roots to detoxify Al in the rhizosphere, thus mak- 
ing R108 plants less tolerance to Al than J. A17 plants. 
In addition, we demonstrated that promoter sequence in 
MtZpt2-l of J. A17 plants contained more response- 
elements than that of R108 plants. Given the regulatory 
roles of MtZpt2-l in response to salt stress, these results 
may account for the greater tolerance of J. A17 plants to 
salt stress than R108 plants. Finally, our results revealed 
that deletion of an YSL gene encoding an iron trans- 
porter in the genome of R 108 plants is likely to impair 
long-distance transport of iron in R108 plants. This re- 
sult may explain the greater sensitivity of R108 plants to 
iron deficiency than J. A17 plants. Taken together, these 
findings demonstrate that analyses of genome variations 
by sequencing can shed important light on differences in 
responses of M. truncatula ecotypes to abiotic stress in 
general and mineral stress in particular. 

Methods 

Plant materials and treatments 

Two Medicago truncatula ecotypes Jemalong A17 and 
R108 were used in this study. Seeds of both ecotypes 
were treated with concentrated sulfuric acid for 8 min, 
and then thoroughly rinsed with water. After chilled at 
4°C for 2 d, seeds were sown on 0.8% agar to germinate 
at 25°C until the radicals were approximately 2 cm. The 
seeds were planted in the same plastic buckets (6 seed- 
lings for both ecotypes per bucket) filled with 2.5 L aer- 
ated nutrient solution. The composition of full-strength 
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nutrient solution is: 2.5 mM KN0 3 , 0.5 mM KH 2 P0 4 , 
0.25 mM CaCl 2 , 1 mM MgS0 4 , 100 uM Fe-Na-EDTA, 
30 uM H3BO3, 5 uM MnS0 4 , 1 \iM ZnS0 4 , 1 |iM 
CuS0 4 and 0.7 uM Na 2 Mo0 4 with pH of 6.0. 

For measurements of the effect of A1C1 3 on root elong- 
ation, 3-d-old seedlings were transferred into solutions 
containing 0.5 mM CaCl 2 with and without 5 uM A1C1 3 
(pH 4.5) for 2 days. Length of primary root was mea- 
sured after treatment with A1C1 3 , and relative root elong- 
ation was calculated. To determine the effect of A1C1 3 
on exudation of citrate from roots, three-week-old seed- 
lings were transferred into solutions containing 0.2 mM 
CaCl 2 with and without 5 uM A1C1 3 (pH 4.5) for 1 days. 
The exudation from the treated roots was collected at 
room temperature without light for 2 hours, and then 
citrate concentration in the exudation solution was de- 
termined by reversed-phase high performance liquid 
chromatography (HPLC) as described previously [56]. 
For measurements of Al content in roots, seedlings of 
the two ecotypes were treated with 5 uM A1C1 3 (pH 4.5) 
for 24 h, and roots were collected for measurement 
of Al. 

Three-week-old seedlings were transferred into solu- 
tions containing 100 mM NaCl or 1 uM Fe-Na-EDTA 
for 5 days. Shoots were collected to measure the content 
of Na + , K + and Fe. 

Measurement of mineral elements 

Plant materials treated with and without mineral stress 
(Al and Na toxicity and Fe deficiency) were harvested 
and dried at 80°C to constant weight. As much as 50 mg 
of dry plant material was weighed and placed in a diges- 
tion tube, and then samples were digested with 6 mL of 
nitric acid and 2 mL of hydrogen peroxide using micro- 
wave system (MARS, CEM). The digest were diluted to 
50 mL. After filtering, the concentrations of Al, Na, K 
and Fe were measured by ICP-AES (Thermo). 

DNA isolation and resequencing 

DNA isolation was carried out using a CTAB (cetyl tri- 
methylammonium bromide) protocol. After quality assay, 
genomic DNA was fragmented randomly. After electro- 
phoresis, DNA fragments of about 500 bp were gel puri- 
fied. Adapter ligation and DNA cluster preparation were 
performed and subjected to 2 x 90 bp paired-end sequen- 
cing on an Illumina Hiseq2000 sequencer. The raw data 
have been submitted to NCBI Sequence Read Archive 
(http://www.ncbi.nlm.nih.gov/sra) and the accession num- 
ber is SRP029924. 

Bioinformatics analysis 

Firstly, adapter contamination in the raw data was re- 
moved. To ensure quality, each base in a read was 
assigned a quality score (Q) by a phred-like algorithm 



[57,58]. The reads which contained more than 50% 
low quality bases (Q < 5) were removed. Using SOAP2 
[59], all reads were aligned with the M. truncatula 
reference genome (Mt 3.5 assembly) [10]. If the ori- 
ginal read could not be aligned onto the reference se- 
quence, the first nucleotides at 5' end and two nucleotides 
at 3' end were deleted, and then aligned onto the ref- 
erence again. If the sequence failed to alignment, two 
more nucleotides at 3' end were deleted. The pro- 
cedure was repeated until alignment was achieved or 
the read was less than 32 bp. The average sequen- 
cing depth and coverage was calculated using the re- 
sults of alignment. 

Structure variations, short indels and SNPs were 
identified by aligning the reads of R108 obtained from 
resequencing to the reference genome of ecotype J. 
A17. In our experiment, the distance of both relevant 
paired-end reads should be about 500 bp. However, if 
the distance and orientation were different from ex- 
pectation after both relevant paired-end reads were 
aligned with the reference genome, the region might 
have variation structures. The types of structure vari- 
ations that can be detected include deletion, insertion, 
duplication, inversion and transposition. SOAPsv was 
used to identify structure variation, and at least three 
paired-end reads were needed to confirm a variation 
structure in the present study. The alignment gaps 
in mapped reads were identified as candidate indels 
using SOAPindel. The maximum gap length was 5 bp, 
and at least three pairs of reads to define an indel. 
On the basis of alignment, polymorphic loci against 
the reference sequence were identified according to 
the following criteria: Q > 20, 3 < Depth < 100 and at 
least 5 bp away from each other. SOAPsnp was used 
in this assay. 

Plant ds-acting regulatory elements were searched by 
the PLACE database [60]. 

RNA isolation and real-time quantitative PCR 

Total RNA was isolated using RNAiso Plus reagent 
(TaKaRa) and treated with RNase-free DNase I (Promega). 
The total RNA was reverse-transcribed into first-strand 
cDNA with PrimeScripf RT reagent Kit (TaKaRa). 

Real-time quantitative PCR (RT-qPCR) was performed 
using ABI Stepone Plus instrument. Gene-specific 
primers of MtAACT (accession No. XM_003627650.1) 
were 5'-GAC ATA GAG AAA GGG ACA-3' and 5'- 
AGG ATA GTA AAT GGG GTT-3'. MtActin (accession 
No. BT141409) and MtGADPH (accession No. XM_ 
003608827.1) were used as internal control with primers: 
(5'-ACG AGC GTT TCA GAT G-3' and 5'-ACC TCC 
GAT CCA GAC A-3') and (5'- AAG GAG GAG TCT 
GAG GGC-3' and 5'-AAC GGC TGC TAG GCT AAT- 
3'). Each reaction contained 5.0 uL of SYBR Green 
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Master Mix reagent (TOYOBO), 0.4 \iL cDNA samples, 
and 0.6 \iL of 10 \xM gene-specific primers in a final 
volume of 10 \iL. The thermal cycle used was 95°C for 
2 min, 40 cycles of 95°C for 30 s, 55°C for 30 s, and 
72°C for 30 s. The relative expression level was calcu- 
lated by the comparative C T method. 
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